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Abstract 

Based on transfer matrix techniques and finite size scaling, we study 
the oriented polymer (self-avoiding walk) with nearest neighbor interac- 
tion. In the repulsive regime, various critical exponents are computed and 
compared with exact values predicted recently. The polymer is also found 
to undergo a spiral transition for sufficiently strong attractive interaction. 
The fractal dimension of the polymer is computed in the repulsive, attrac- 
tive regimes and at the spiral transition point. The later is found to be 
different from that at the collapse transition of ordinary self-avoiding walk. 
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1 Introduction 



A polymer chain describes a self avoiding walk (SAW). In the high temperature 
phase, one generally ignores the Van Der Waals attraction between monomers of 
the polymer chain and study the polymer problem as the n — )• limit of the 0{n) 
model[l]. In 2d, this has been the basis for calculation of many critical exponents 
that have important physical meaning[l, 2, 3, 5], and also starting point for 
many interesting extensions of the polymer problem such as polymer in the low 
temperature phase where monomers' attractions are taken into account [4, 5, 6, 
7]. However, much less attention has been paid to the study of polymer chain 
with orientation, ie. polymer with an 'arrow' that runs along the chain[8, 9]. 
Physically, one can imagine such orientation as arising from dipole moments of the 
monomers that make up the chain. Clearly, if we only consider the interactions 
between monomers to be the excluded volume ones, then the orientation does 
not play any role and there is no difference between the oriented and ordinary 
polymers. Indeed, an oriented polymer chain can be considered as the n — )• 
limit of the complex 0{n) model whose partition function is given by, 

Y.I[i^ + ^s,-s* + xs:-s,) (I.I) 

{5} <«J> 

where Si is a complex n-component vector at site i with * denotes complex 
conjugation, <> denotes nearest neighbor lattice sites, and x is the fugacity of 
the monomer. The complex 0{n) model is the same as the 0{2n) model, hence 
the oriented SAW is in the same universality class as the ordinary SAW. 

An interesting extension of this model consists in introducing interaction be- 
tween the monomers that depends on their relative orientation[9]. More precisely, 
consider two monomers that are in close encounter, (for a square lattice, for ex- 
ample, this can be defined by neighboring lattice edges being occupied by the 
monomers) their relative orientation can either be parallel or anti-parallel. We 
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assign different Boitzmann weights f or e ° to these two types of ciose encounters 
(see fig. f). Note that Ao here differs from that in [9] by a sign. 

Such an orientation dependent interaction has been identified in the context 
of the complex 0{n) model as a current-current interaction[9]. The complex 0{n) 
model possesses a U{1) symmetry associated with S — > e^"S, the corresponding 
current = ^.S** • df^S — S ■ d^S*) can be identified as the tangent vector along 
the oriented chain[8]. Perturbing the complex 0{n) model by a current-current 
interaction, A / J Jdzdz, corresponds to introducing an orientation dependent in- 
teraction for the monomers of the polymer chain. 

This model includes the ordinary SAW as a special case where Ao = 0. For 
Ao ^ 0, parallel close encounters are suppressed. While for Ao ^ parallel 
close encounters are preferred, the polymer has the tendency to coil up to form 
a spiral. We shall denote these two phases respectively as the repulsive and 
attractive regimes. For some intermediate Ao > we expect a transition to take 
place. This transition, denoted as the spiral transition, should be different from 
the usual collapse transition (or the 0{0 ) point) of the ordinary SAW[5, 6, 7] as 
one can readily see that in this case it is the self-trapping configurations that are 
responsible for the collapse. We therefore expect a new multi-critical point from 
this spiral transition. 

On the other hand, in the repulsive regime, regarding the current-current in- 
teraction as a perturbation from the ordinary SAW, which belongs to this regime, 
and using the fact that the perturbation is truly marginal[10], critical exponents 
can be determined exactly as functions of A (the renormalized Ao)[9]. For exam- 
ple, the operator which corresponds to the source of a L-\eg oriented SAW with 
U{1) charge Z, (ie. L arrows originate from the source) scales with exponent 



xl{\) = xl{0) - 2tt\L^ 




where 



9L^ 
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is the critical exponent for the ordinary SAW. Also in this case, the fractal di- 
mension of the SAW is still given by that of the ordinary SAW, Dp = |, since 
the energy operator Si ■ S* + S* ■ Sj is associated with self avoiding ring, which 
can not have parallel close encounter. 

In this paper, we present a numerical study of this interacting oriented poly- 
mer problem based on transfer matrix and finite size scaling techniques. We 
mainly consider oriented SAWs on a square lattice with periodic boundary con- 
dition. For the repulsive regime, we verify the various critical exponents xl[X) 
conjectured, and demonstrate the presence of a spiral transition as one increases 
Xo- On the spiral transition point, exponents and 7 are determined numerically. 

The basic quantities that we are interested are the connectivity constant and 
the critical exponents. They can be determined from the spin-spin correlation 
function. The spin-spin correlation function in the high temperature expansion 
can be expressed as a sum over all configurations of oriented SAW with source 
and sink at the locations and R of the two spins respectively, such that each 
configuration is weighted by the length and the number of parallel close encounters 
as given below: 

< 5^(0) • S*{R) >= J2 x'(^)e^°"(^) . (1.4) 
c 

Here l[C) and m[C) are respectively the length ( number of steps ) and the number 
of parallel close encounters of the oriented SAW C. Note that in this model we 
take the weight of anti-parallel close encounter to be 1 as this choice ensures that 
the critical fugacity is independent of Xo in the repulsive regime[9], to which we 
shall mainly restrict our attention. 
The expansion can be rewritten as 

< 5^(0) • S*iR) >= ZiiO, R)x' (1.5) 

where 

Z,(0,i?) = ^cu,,„(0,i?)e^°" (1.6) 
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is defined as tlie partition function of a polymer cliain. Hence tlie correlation 
function is a generating function of the polymer partition function. In the vicinity 
of the critical point, x — )• x^", 

< 5^(0) • S*{R) >~ ^^e-l^l/« \R\ ^ oo (1.7) 

where Xi is the exponent given in eqn.( 1.2) and the correlation length ^ diverges 
^1 ^ . The constants x^, Xi and v have geometrical meanings as can be 
seen by inverting the above generating function. One readily shows that the SAW 
partition function has the asymptotic form 

Z,(0, R) ~ (x=)-'r(^)-^-2''F(|i?|/r) /, \R\ oo (1.8) 

where 

7(A) = //(2-2xi(A)) (1.9) 

and F is some scaling function. Hence, x'^~^ is the connectivity constant that 
depends on the type of lattice, which by our choice of the interaction energy, 
does not depend on A (Ao). For a square lattice, the value has been determined 
with high precision to be 0.3790528(25)[11]. The thermal exponent which is also 
independent of A is equal to the inverse of the fractal dimension of the model 
since the average radius of gyration 

< \R\ >= J2 \R\Zi{o, R)l E ^'(0' R) (1-10) 

R R 

behaves asymptotically as 

< 1^1 l^oo. (1.11) 

So geometrically, (x^)"^ and 7 determine the number configurations of a polymer 
with given length, while v determines the average size. 

Besides the spin-spin correlation function, one can also consider correlation 
function of an operator cf)^ with U{1) charge L that corresponds to insertion of the 
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source of a L-leg oriented SAW. In this case, the correlation function generates 
the 'watermelon' diagram consisting of L mutually and self avoiding oriented 
walks tied together at their extremities. 

2 Polymer on a strip 

Consider a polymer on a strip of finite width n. The generating function of 
Z;(0, i?), and hence the constants x^, 7, and can be obtained from the transfer 
matrix. The method employed is similar to that used in the study of ordinary 
polymer [11, 3, 5]. A state is still characterized by monomers' configurations on 
a column of the strip as well as the connectivity to the left of these monomers. 
But in our case, monomers and their connectivity carry orientations. As a result, 
besides the usual non-crossing constraint, connectivity has to satisfy additional 
constraint which may be regarded as the non-crossing constraint for connectivity 
to the right of the monomers. As an example, the column of monomers and their 
connectivity shown in fig. 2 can not represents section of a single SAW. 

To keep track of the monomers' interactions, one has to consider monomers' 
configurations on a column of horizontal and vertical edges, this together with the 
fact that connectivity carries orientation vastly increases the number of possible 
states. In most cases we manage to diagonalize transfer matrix of dimension of 
the order of 7000 X 7000, which for the one polymer sector is equivalent to a strip 
of size 7 with translational symmetry taken in account. We list some of their 
dimensions in table 1 where SL,n denotes dimension of the translational invariant 
sector with L polymers and n is the width of the strip. 

For large the generating function eqn.( 1.4) is given by the largest eigen- 
value A of the transfer matrix of the one polymer sector as A'^' and this defines 
a correlation length Ao) by 
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The critical point is defined as the smallest where the correlation length di- 
verges (ie. A(a;^,Ao) = 1). Keeping only the dominating term in eqn.( 1.8), the 
free energy is simply given by 

/4A„) = -log<(A„) (2.2) 

from which one can determine all thermodynamics quantities. In particular, the 
susceptibility with respect to A can be obtained as 

i! 

dXt 



Cn{\o) = -^,\0gxl (2.3) 



We present in fig. 3 a plot of this susceptibility against Ao for various strip size 
n ranging from 3 to 7. This susceptibility exhibits a peak near Ao ~ 1, signaling 
a possible phase transition. 

Using the partition function eqn.( 1.8), the average radius of gyration of a 
polymer chain of length / can be expressed as 

Since there is only one polymer in the area n < \Ri\ >, the average density of the 
chain is thus given by 

1 d 

PniK) = -w] log A«, Xo) . (2.5) 

n o log X 

This density is plotted in fig. 4, it increases with Ao, where Ao ~ 1 divides the 
low and high density regimes. Since polymer in the repulsive (attractive) regime 
has low (high) density, Ao ~ 1 is roughly the spiral transition point. This change 
of density can also be seen from the coefficient of expansion defined as 

tn{X.) = (2.6) 

In fig. 5, We again notice a peak around Ao ~ 1. 

From these figures, we see strong evidence of a singularity near Ao ~ 1 building 
up as n — )• oo. In table 2 locations of successive maximums of c„ and t„ are listed 
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and these values converge to 1.1 7 it 0.2 and 1.12 it 0.2 respectively. It is very likely 
that this is the location of the spiral transition point. We shall denote this point 
as A*. 

So far we have concentrated on the critical line, but as in the case of ordinary 
polymer, one can give physical meaning to the entire (x, Ao) plane by considering 
the polymer as a chain with flexible length controlled by the fugacity x, the 
analysis is essentially similar to that of the ordinary SAW[3, 5]. In summary, 
every point in the (x, Ao) plane is accessible to the polymer chain, the critical line 
X = x'^[Xo) now corresponds to the isobar where the pressure vanishes. Polymer 
in the region below (above) this line is subjected to negative (positive) pressure 
and is swollen (dense). It can also be shown that crossing this critical line in the 
attractive regime one encounters a second order phase transition, while crossing it 
in the attractive regime one encounters a first order phase transition characterized 
by a jump in the density defined as 

1 d 

Pn{x, Xo) = -TT, log A(x, Ao) (2.7) 

n o log X 

in this grand canonical picture. 

3 Thermal exponent u 

The correlation function (^(x, Ao) diverges at the critical line x = x'^ . For given Ao, 
the divergence is characterized by the exponent v. In the repulsive regime, the 
thermal exponent v is expected to be equal to that of the ordinary SAW ie. 3/4. 
While in the attractive regime, because the polymer is dense, it is expected to be 
equal to 1/2. Note that this value is also expected from the first order transition 
nature of the critical line in this regime[12]. The spiral transition point, which 
divides these two regimes, has v = v^^ assumes some intermediate value. Near x^, 
the correlation length defined in eqn.( 2.1) has the following finite size scaling 
form 

in{x,\o) ^ nF^{{x - x^)n^l'') n > 1, x - < 1. (3.1) 
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From eqn.( 2.5) and log A = — it is clear that 

^^^2) (3.2) 



'log(/^n+i«+i, A)//)^«, A)) 



y log(n + 1 jn) 

converges to ;/ as n — )• oo. 

In fig. 6 we plot successive values of Vn against Xo- The various lines for 
different strip sizes intersect almost at the same point at Ao ~ f where spiral 
transition takes place. Below this point, the thermal exponent converges to 3/4 
rather slowly, while above this point, v decreases to 1/2 as expected from earlier 
prediction. 

An alternative method to determine this exponent uses the phenomenological 
renormalization relation 

-iniXnAo) = —im{x'l,,\o) n, m > I (3.3) 
n m 

which is again based on the finite size scaling assumption eqn.( 3.1). Note that 
the above relation is also used in the first order transition line where Ao > A*, this 
is justified by the fact that the correlation length diverges at the critical point. 
The extrapolated critical value for large n obtained from this relation can 
be checked to be in agreement with that obtained from the jump of the density 

Pn{x, Ao). 

In the numerical calculation we take m = n + 1 in eqn.( 3.3) where convergence 
is fastest to determine successive critical line x = x^. The thermal exponent can 
be determined from the derivative of the correlation length by extrapolating Vn 
defined now as 

The result shown in fig. 7 again shares the same features of the previous figure: 
Thermal exponents v extrapolated in the attractive and repulsive regimes are 
compatible with that obtained before but in this case the convergence is slightly 
faster. Furthermore, various lines cross also at Ao ~ 1. 
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Successive obtained from eqn.( 3.3) are shown in fig. 8. Indeed, the ex- 
trapolated for large n is not sensitive to A in the section of the repulsive 
regime where Ao < A* as expected from the choice of the relative weight for the 
monomers' interactions. So indeed the critical fugacity is equal to that of the 
ordinary SAW 

x=(A„) = 0.3790528(25) . (3.5) 

Another interesting feature of these thermal exponents is the crossings of the 
various lines, which converge to the spiral transition point as n — )• oo. The 
corresponding thermal exponents and Ao values for various crossing of the lines 
are listed in table 3. Due to the small number of data and system sizes we have, 
it is difficult to obtain any reliable extrapolated values. Nonetheless, from both 
sets of data, it is likely that 

a: ^ 1.1 ±0.2 

~ 0.63 ±0.05 ^ ^ 

It should be remarked that since the free energy (eqn.( 2.2)) is always greater 
than the free energy associated with a tightly wound spiral (which is given by 
Ao), and that in the entire repulsive regime up to the spiral transition point is 
given by eqn.( 3.5) independent of Ao. We have 

log{xT' > Ao for Ao < a: . (3.7) 

This implies the following bound 

A* < 0.970080(7) . (3.8) 

The values for A* deduced from tables 2, 3 are very close to this bound. The 
fact that they exceed this bound can be attributed to the small system sizes we 
can handle and hence the unreliability of the extraploted values. Note that the 
estimated A* is very close to the value of the bound ie. log(a;^)~^. It is tempting 
to conjecture that indeed the A* is given by the bound ie. the value obtained from 
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the critical lugacity ol ordinary SAW . II this is indeed the case then this implies 
that there is no 0(1) excitation from the tightly wound spiral, which suggests 
that the spiral transition is first order and quite different from the 6 point for the 
collapse of ordinary SAW. We hope to pursue this issue in the future. 

4 Conformal Weights 

In the repulsive regime the exponent xl[X) that characterizes the scaling behav- 
iors of the iv-leg oriented SAW can be related to that of the ordinary SAW. The 
latter has been identified as conformal weight of a twisted N = 2 supersymmetry 
theory with k = 1[14]. Thus, the oriented SAW in this regime is described by a 
marginal perturbation of the twisted N = 2 theory in the continuum limit [15]. 
This exponent can likewise be computed using the transfer matrix technique. 
As in eqn.( 1.4), we consider the correlation function of the L-leg operator (pL 
on a finite strip of width n and define a correlation length ^L^n from the largest 
eigenvalue of the corresponding transfer matrix as 

AL,n = e-'/^^'" . (4.1) 
From conformal field invariance argument[f 7], XL^n{^o) defined as 

converges to xl[X) as n — )• oo. Here the correlation length is computed at the 
critical point defined using relation similar to ( 3.3). 

In tables 4,5,6,7, we list the conformal weights computed at Ao = —3, —2, —1,0 
for various L and strip widths n. The extrapolated values for n — )• oo are also 
given in the same tables. It is however not straight forward to compare these 
results with eqn.( 1.2) since the relation between A and Ao is not known except 
at the point where they both vanish. For this ordinary SAW point, the data are 
listed in table 7, here we clearly see that the extrapolated values agree very well 
with the exact values. In general convergence is better for sectors with small L. 
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In fig. 9, we plot successive values of Xi^n{^o) against Ao for L = 1. Certain 
features are obvious from the figure: First, the conformal weight increases as Ao 
decreases as predicted by eqn.( 1.2). Note that A and Ao should have the same 
sign. Second, the crossings and minimums (for large n) of the lines seems to 
converge to the point A* given in eqn.( 3.8). Furthermore, the value of a;i^„(A*) 
at the crossing (minimum) converges to zero as n increases. This supports the 
conjecture that the spiral transition takes place at Xi[X) = 0[9]. From eqn.( 1.2), 
we deduce the renormalized A* to be 5/967r. 

The vanishing of a;i(A*) also implies that the exponent 7 at the spiral transi- 
tion point is given by 

7 = 2z/^ ~ 1.26±0.1 (4.3) 

using eqn.( 1.9). It is intriguing that such a relation between 7 and also exists 
for the ordinary SAW at the 6(6 ) point [13]. It should be stressed that 7 obtained 
in this way is certainly not rigorous, it would be desirable to check it against other 
methods. For example, the computation of the generating function 

G{x, Ao) = ^ < 5^(0) • s{R) >~ \x - (4.4) 

approaches x^. 

One could also compute quantities that are independent of A. For example, 
the following ratio 

should converge to 

MA)+izi2^/xy 

x^'(A) + l/12 \L'J ^ ■ ' 

which depends only on the ratio of the number of 'legs'. However, because the 
size of our strip is small, successive r[-"|;'s do not form a monotonic sequence and 
extrapolation is difficult. We therefore compute the ratio 2.' directly from the 
extrapolated values of the conformal weights given in tables 4,5,6,7. The results 
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are presented in table 8 below. Convergence is better near Ao = 0, nevertheless, 
we see evidence that the data agree to the exact values which are independent of 
A.. 

5 Conclusion 

In this numerical study, we verify several predictions and conjectures made in 
[9] regarding the repulsive regime. In particular, the critical exponent xl[X) is 
checked and thermal exponent p is found to have the value 3/4 independent of 
A. Much less is known about the attractive regime except that p is found to be 
1/2 consistent with the fact that the polymer in this regime is dense. 

The numerical data also give strong evidence that a spiral transition occurs 
as the current-current coupling Ao becomes large. The location of this transition 
is determined with less accuracy to be 



which is very close to the bound given in eqn.( 3.8). This could shed light on 
the nature of the spiral transition. On this transition point, we also deduce the 
exponents 



which are different from that of the usual collapse transition of the ordinary SAW. 
This spiral transition certainly deserved further investigation. 

Besides the critical line x = x^(A), the region x > x'^[X) also deserves fur- 
ther investigation. For the ordinary SAW, it has been known for a long time 
that this region is critical and described by the low temperature phase of the 
(9(0) model[16], exact expressions for the critical exponents that have geometri- 
cal meaning have also been given. It would be interesting to study the effect of 
the current-current perturbation of this phase. We hope to pursue these questions 
in the future. 



A* ~ 1.1 ±0.2, 




= 0.63 ±0.05 
7 = 1.26 ±0.10 
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Table Captions 

• Table 1 Dimensions ol various polymers sectors. 

• Table 2 A* obtained from maximums of successive c„ and t„ respectively. 

• Table 3 A* and obtained from crossings of lines in fig. 6 and fig. 7 respec- 
tively. 

• Table 4 Conformal weights at Ao = —3 from different strip widths and L. 

• Table 5 Conformal weights at Ao = —2 from different strip widths and L. 

• Table 6 Conformal weights at Ao = — 1 from different strip widths and L. 

• Table 7 Conformal weights at Ao = from different strip widths and L. 

• Table 8 Ratio of conformal weights at various Ao's. 
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Figure 2: Monomers and Connectivity. 
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Figure 4: Density versus Ao- 



20 




21 




-2-10 12 3 Ao 



Figure 6: Exponent v against Ao for different strip widths n. 
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n 




S2,n 


S3,n 


S4,n 


S5,n 


S6,n 


2 


3 


1 










3 


11 


3 


1 








4 


45 


17 


3 


1 






5 


195 


78 


20 


3 


1 




6 


881 


402 


122 


26 


3 


1 


7 


4121 


2017 


689 


171 


29 


3 


8 


19831 


10470 


3877 


1114 


231 


35 


9 






21630 


6802 


1656 


302 


10 










11140 


2381 



Table 1: Dimensions of various polymers sectors. 



n 


3 


4 


5 


6 


7 


extrapolation 


K 


1.85 


1.75 


1.57 


1.41 


1.32 


1.17(20) 


K 


1.82 


1.65 


1.53 


1.41 


1.33 


1.12(20) 



Table 2: A* obtained from maximums of successive c„ and t„ respectively. 
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n 


3,4 


4,5 


5,6 


K 


1.166 


1.128 


1.111 




0.628 


0.642 


0.641 


K 


1.016 


1.084 


1.087 




0.645 


0.627 


0.622 



Table 3: A* and Vx obtained from crossings of lines in fig. 6 and fig. 7 respectively. 



n 










2 


0.12918 


2.52743 






3 


0.12553 


2.08935 


5.68672 




4 


0.12275 


1.05229 


5.31323 


10.01097 


5 


0.12022 


0.92181 


4.12591 


9.80083 


6 


0.11810 


0.85923 


2.61618 


8.58174 


7 


0.11640 


0.82220 


2.32699 


6.78043 


8 








4.77810 


9 








4.34188 


extrapolation 


0.116(8) 


0.736(10) 


1.731(9) 


3.91(11) 



Table 4: Conformal weights at Ao = —3 from different strip widths and L. 
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n 


^^1,71 




^^3,71 




2 


0.12918 


1.89081 






3 


0.12479 


1.61188 


4.25433 




4 


0.12186 


1.01428 


4.03999 


7.56325 


5 


0.11936 


0.90026 


3.31607 


7.41350 


6 


0.11732 


0.84433 


2.50402 


6.66471 


7 


0.11570 


0.81079 


2.25849 


5.61332 


8 








4.57459 


9 








4.20151 


extrapolation 


0.114(9) 


0.724(8) 


1.720(9) 


3.33(11) 



Table 5: Conformal weights at = —2 from different strip widths and L. 



^^1,71 ^^2,71 ^^3,71 ^^4,71 



2 


0.12918 


1.25419 






3 


0.12272 


1.13442 


2.82193 




4 


0.11926 


0.91246 


2.78875 


5.01677 


5 


0.11682 


0.83990 


2.47811 


5.02618 


6 


0.11499 


0.80148 


2.20838 


4.73237 


7 


0.11359 


0.77747 


2.06935 


4.34551 


8 








4.01605 


9 








3.81281 


extrapolation 


0.112(9) 


0.703(8) 


1.690(8) 


3.22(9) 



Table 6: Conformal weights at Ao = — 1 from different strip widths and L. 
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n 




^2 77, 


X3n 




Xs 


2 


0.12918 


0.61757 








3 


0.11651 


0.65695 


1.38954 






4 


0.11076 


0.66641 


1.49351 


2.47029 




5 


0.10809 


0.66892 


1.54045 


2.63885 


3.85982 


6 


0.10676 


0.66948 


1.56454 


2.73158 


4.09298 


7 


0.10601 


0.66946 


1.57813 


2.78720 


4.23650 


8 








2.82302 


4.33011 


9 








2.84573 


4.39399 


extrapolation 


0.1041(6) 


0.669(7) 


1.603(7) 


2.92(8) 


4.598(7) 


exact 


0.10417 


0.66667 


1.60417 


2.91667 


4.60417 



Table 7: Conformal weights at Ao = from different strip widths and L. 





-3 


-2 


-1 





exact 


r2,i 


4.11(14) 


4.08(13) 


4.02(16) 


4.01(12) 


4.00 


rs,! 


9.10(27) 


9.12(27) 


9.08(31) 


9.00(21) 


9.00 


r-A,! 


17.44(89) 


17.28(84) 


16.90(83) 


16.02(66) 


16.00 


r3,2 


2.21(6) 


2.23(5) 


2.26(5) 


2.24(6) 


2.25 


r4,2 


4.24(15) 


4.22(13) 


4,20(12) 


3.99(12) 


4.00 


r4,3 


1.91(7) 


1.89(8) 


1.86(8) 


1.78(9) 


1.78 



Table 8: Ratio of conformal weights at various Ao's. 
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